Hydrodynamic theory of transport in doped graphene 
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We study non-linear dc transport in graphene using a hydrodynamic approach and conclude that 
in clean samples the drift velocity saturates at a weakly density-dependent value Vsat ~ 10^ cm/s. 
We show that saturation results from the interactions between graphene's Dirac quasi-particles and 
both acoustic and optical phonons. Saturation is accompanied by substantial electron heating and 
is not reached at realistic driving fields in moderately or strongly disordered samples. We find that 
it is essential to account for interactions among graphene's Dirac quasi-particles, which increase the 
linear response resistivity at high temperatures or low densities. 

PACS numbers: 71.35.-y,73.21.-b,73.22.Gk,71.10.-w 



I. INTRODUCTION 

The unique electronic and thermal properties of 
graphene two-dimensional electron systems make them 
promising as potential building blocks for future elec- 
tronic devices.-'^ The gapless Dirac-like spectrum of 
graphene presents an obstacle to logic device applica- 
tions, but is also partly responsible for large quasiparticle 
velocities which are advantageous in analog and radio- 
frequency devices.— The feasibility of these applications 
is dependent mainly on graphene's non-linear electrical 
response properties which we address in this paper using 
a hydrodynamic approach. 

The response of graphene to an external electric field 
is determined by the interaction of its Dirac quasi- 
particles with impurities, with phonons, and with each 
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For samples on substrates, the resis- 



tivity contribution due to scattering from the phonon 
modes of graphene is unimportant in the linear response 
regime even at room temperature. The linear resistivity 
is generally believed to be limited by elastic scattering off 
Coulomb impuritiesi^ii^'i^ Even in annealed suspended 
graphene sheets, which have dramatically weakened elas- 
tic scattering, the phonon-limited resistance is irrelevant 
in currently available samples because it is small com- 
pared to the quantum resistance. ^^'^^ At high-fields, how- 
ever, electron-phonon (e-ph) interactions are essential for 
two reasons. First, the phonon modes are the only dis- 
sipative channels through which the electrons can lose 
the energy that they acquire by flowing in the presence 
of a high electric field. Second, the electronic tempera- 
ture rises and the drift velocity u increases with the elec- 
tric field; both changes enhance e-ph scattering, making 
it much more efficient for momentum relaxation. Ex- 
trinsic phonons localized near the substrate surface can 
play an important role even in the linear regime.^'-'^^ How- 
ever the significance of these phonons is sensitive to the 
specific experimental system, in particular to the type 
of substrate and its distance from the graphene sheet. 
We therefore do not account explicitly for these extrinsic 
phonon modes and instead limit ourselves to an explana- 
tion of how they may be straightforwardly incorporated 



when the parameters appropriate for a particular exper- 
imental system are known. 

One very important consequence of the high elec- 
tronic temperature in the non linear regime is enhanced 
electron-electron (e-e) scattering. As recent spectroscopy 
measurementsiiii^ demonstrate, e-e scattering is the 
dominant scattering mechanism at high electronic tem- 
peratures. High temperatures are inevitable in the non 
linear regime because of Joule heating. In the linear 
regime e-e interactions dominate either when the sample 
is hot or when it is disorder-free. When it is dominant, 
the influence of e-e scattering on transport can not be 
considered perturbatively. In this paper we use a theo- 
retical approach to non-linear dc transport which exploits 
rapid e-e collisions by using a hydrodynamic theory. One 
important advantage of this theory is its simplicity and 
physical transparency. As we explain below, in a hydro- 
dynamic theory the non-equilibrium system is character- 
ized by only three parameters: the chemical potential fi, 
the electronic temperature and the drift velocity u. 

Our paper is organized as follows. In section we de- 
rive the hydrodynamic equations. Using these equations 
we study the linear dc transport in graphene in section 
mil We then consider the non- linear case in section [TVl 
We first focus on a clean system and then discuss the role 
of disorder. Finally we summarize our findings in section 
IVl As we explain below, the hydrodynamic description 
breaks down in the neutral regime when \fi\/Te <ti 1. We 
therefore restrict our study to doped systems. Because 
the Dirac model for graphene is perfectly particle-hole 
symmetric we can restrict our attention to electron-doped 
systems without loss of generality. 



II. HYDRODYNAMIC THEORY FOR DOPED 
GRAPHENE 

The Boltzmann theory provides a simple but faithful 
description of transport in many electronic systems. In 
graphene, the validity of this semiclassical transport the- 
ory is well established for the doped systems we study. 
The backbone of Boltzmann transport theory is the dis- 
tribution function /ka(r,t) defined as the occupation 
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probability of the Bloch state in band a with crystal mo- 
menta k at position r and time t. All physical quantities 
can be expressed in terms of /. The distribution function 
is determined by requiring that it satisfies the Boltzmann 
equation: 



{dt + Vk„ • V + eE • Vk) /ka(r, t) = SeL + Se 



(1) 



where Vka is the band velocity, E is the electric field, SeL 
is the electron-lattice collision integral which accounts for 
electron scattering by phonons and disorder, and See is 
the e-e collision integral which accounts for e-e scattering. 

Further simplification of the theory is possible in the 
hydrodynamic regime when the e-e scattering time Tee 
is considerably shorter than the e-ph scattering time Tph 
and the impurity scattering time r^, i.e. when 



'^ee '^ph 5 '^i 



(2) 



This separation of time scales implies that, to leading 
order in both Tee/n and Tee/Tph, rapid e-e collisions are 
able to establish a drifting Fermi distribution function^^ 



exp 



Cka - u(r,t) ■ k - /i(r,t) 
Te(r,t) 



(3) 

where eka is the energy dispersion of band a. This form 
of distribution function satisfies See{f) — 0. The drift 
velocity u can be non-zero because e-e scattering does 
not relax momentum. The existence of the three hydro- 
dynamic functions Tg , u and /i is a direct consequence 
of the conservation of energy, momentum, and particle- 
number in e-e collisions. 

Specializing to dc transport in graphene, the Boltz- 
mann equation reduces to 



eE • Vk/,^„ = SeL + Se 



(4) 



In Eq. Q the hydrodynamic functions are independent of 
time and position, and eka = avk with v being the band 
velocity of graphene. In what follows we interchangeably 
use a = c,v and a = to label the conduction and 
valence bands of graphene. 

The three hydrodynamic parameters give a full de- 
scription of the non-equilibrium state. Their values are 
fixed by the aforementioned conservation laws. Multi- 
plying Eq. ^ by k, summing over it and the band index, 
and using X^ka = implied by the conservation of 
momentum in an e-e scattering event we obtain the force 
balance equation: 



enE = V 



where 



V^-gY^kSeL ifU- 



(5) 



(6) 



ka 



Here 5 = 4 accounts for the spin and valley degeneracies 
in graphene. The equality in Eq.® expresses the steady 



state balance between the momentum acquired by the 
charge carriers due to the electric field and the momen- 
tum lost by scattering off phonons and impurities. 

Similarly by multiplying Eq.® by eka, summing 
over momenta and over the band index and using 
X^ka ^kaSee = wc obtain the energy balance equation: 



enE • u = Q 



where 
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E 



£kaSeL (/ka) ■ 



(7) 



(8) 



Q is positive when Te > T^. In dc transport the energy 
gained by carriers due to drift in an electric field must be 
balanced by energy lost to the phonon bath. 

Finally, the third equation necessary to fix the values 
of the hydrodynamic variables 



5L[/kc-(l-/kJ 

k 



(9) 



follows from number conservation. Eqs. (j5l7p and ^ de- 
termine the hydrodynamic parameters /i, Te and u , given 
the values for the electric field E, the lattice temperature 
Tl and the density n. 

In subsequent sections we use the hydrodynamic equa- 
tions to study dc transport in graphene. We start in the 
next section by considering the linear electrical response 
of graphene, before turning in the following section to 
the full non-linear response. 



III. LINEAR RESPONSE 

When high currents are driven through a graphene 
sheet it is heated i.e. > Tl. The chemical poten- 
tial is then reduced relative to its equilibrium value to 
maintain a fixed electronic density. Symmetry consider- 
ations imply that inversion of the electric field inverts the 
drift velocity, but renders Te and ^ unchanged. There- 
fore in the linear response regime, i.e. to first order in E, 
Te and fi retain their equilibrium values and the energy 
balance equation ([7]) and the number equation ^ are 
satisfied identically. The drift velocity u follows from the 
momentum balance equation ([5]). 

At physically relevant temperatures, momentum loss 
in graphene is thought to be primarily due to long range 
Coulomb scatterers and secondarily due to interactions 
of electrons with longitudinal acoustic phonons. We 
find, in accord with experiment, that the momentum 
loss rate due to the energetic intrinsic optical phonons 
is negligible. Above the Bloch-Griineisen temperature 
Tbg ~ 2cv'7m scattering by acoustic phonons is quasi- 
elastic due to the large mismatch between the sound ve- 
locity c and graphene's band velocity v. For elastic scat- 
tering 



P = .g5]k(/k"„-/«JVF, 



k,p 



(10) 
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where Wk.p is the transition rate between states k and 
p. Expanding to hnear order in u and using 

(11) 



cos 6p Wk,p — cos 9] 

p 



p 



where 6* = 6'k — 6'p is the relative angle between the in- 
coming and outgoing momenta, we find that 



enE : 



u 

-^2 



2tT TkVk 



\-d. (f^ - f^) (12) 



(0) 



from which the drift velocity u readily follows. Here fj^^ 
is the equilibrium Fermi distribution function, Vk = dk^k 
is the band velocity and r^T^ = X]p(l ~ cos0)Wk,p is the 
elastic (transport) scattering rate. The resistivity 
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follows from Eq. and from the expression for the cur- 
rent 



I = e2^Vka/ka = eun. 

ka 



(14) 



To illustrate the influence of e-e interactions on the 
resistivity we compare pee to 



Po = 



(15) 



the resistivity obtained directly from the Boltzmann 
equation when e-e interactions are neglected. At zero 
temperature the hole density in the valence band vanishes 
and the resistivity is not modified by e-e interactions. At 
finite temperatures, however, pee is always larger than 
Po- Note that the resistivity expressions, (|13p and p^ . 
assume only quasi-elastic scattering and isotropy and are 
valid irrespective of the energy dispersion. 

It is instructive to express the resistivity in graphene as 
the sum of the residual resistivity p^-^^ and the acoustic 
phonon induced resistivity p^P'^\ The residual resistiv- 



jk/uQ for 



the momentum relaxation time associated with Coulomb 
scatterers. Here Ug = (yre^/e) where rii is the im- 
purity concentration and e is the dielectric function.^ It 
follows from Eqs. (|13ll5p that 



ity follows from Eqs. (|13ll5p by setting 



4i) 



Po 



nh 



(16) 



where pi*'' = Uq/ [e^w^(ne -|- rih)] . The dependence of p^'^ 
on the number of electrons rie and the number of holes 
nil is qualitatively changed due to e-e interactions. As 
the temperatures is raised both n,, and Uh increase while 
the total density n — — Uh remains fixed. Thus e-e 
interactions change p*-*' from being a monotonic decreas- 
ing function of temperature to a monotonic increasing 



function of temperature. The ratio pie /Po^ plotted in 
figure [TJ 

As the temperature is increased the influence of 
phonon scattering on the resistivity becomes more im- 
portant. The momentum relaxation time associated with 

acoustic phonons is t'"^^^ = v/kC(Ti^) where C{T) = 
'D^T/2pc^ , and V is the deformation potential^. It there- 
fore follows from Eqs. p^ and that above the Bloch- 
Gritneisen temperature 



p(f)=pr'H(M/T). 



(17) 



Here po^''^ = ■kC{T) / e^v^ is the phonon induced resistiv- 
ity in the absence of e-e interactions and 

x^dx [h{x, z) ~\- h{x, —z)] 



n{z) 



(18) 



where h{x,z) = [exp(a:: — z) + 1]^^. The function tL 
is plotted in figure [1] As expected from our previous 
discussion Ti approaches unity in the p/T — » oo limit. In 
the opposite limit Ti w 5.91 (T/p) . 




FIG. 1: Increase in resistivity due to e-e interactions. 

Why is it that strong e-e interactions increase the 
resistivity? When e-e interactions are neglected elec- 
trons and holes contribute additively to the current. For 
nearly neutral systems, this Boltzmann theory property 
explains the increase of the conductivity as the temper- 
ature is raisedi^ii^. In the opposite hydrodynamic limit, 
strong e-e interactions enforce a common drift velocity 
for all momenta and for both valence and conduction 
bands. The currents in the two bands then flow in oppo- 
site directions, resulting in a smaller net current or equiv- 
alently in a higher resistivity. The diff'erence between the 
magnitudes of the two counter-flowing currents decreases 
as r/|p| is increased. These counterflow currents are 
reminiscent of the Coulomb drag effect in electron-hole 
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bilayersSS. when the two layers are contacted simultane- 
ously. The hydrodynamic description breaks down for 
an electron doped system once the momenta acquired by 
the valence band holes due to the electric field exceeds 
the e-e induced momenta transfer between the bands. 
Therefore our theory is invalid in the neutral regime when 

We now turn to study the non linear dc electric re- 
sponse of graphene. 

IV. NON LINEAR RESPONSE 

Strong electric fields drive the system out of equilib- 
rium. Far from equilibrium the chemical potential is re- 
duced relative to its equilibrium value and the electronic 
temperature is higher than the lattice temperature. To 
find the hydrodynamic parameters which characterize the 
non-equilibrium state, we solve the three coupled hydro- 
dynamic equations (|5l7p and © . We start by simplifying 
the general expressions (jH) and ([8]) for the momenta loss 
V and energy loss Q. In the following we assume that 
E II X. 

A. Momentum loss 

The momentum loss rate of the electronic system 

kct P7 

(19) 

has contributions due to scattering by disorder and by 
acoustic and optical phonons. To calculate V we must 
evaluate the transition rate W and the momentum loss 
rate for each scattering mechanism. 

For elastic collisions expression can be simplified: 



kpa 



(20) 



Tk 



■cos9k [/k"c - (1 - /kJ] 



Obviously, elastic scattering can not induce inter-band 
scattering in graphene . The momentum loss rate due to 
scattering off Coulomb impurities. 



w(l-/32)3/2 Z-^ 



(21) 



is obtained by substituting r^*-* in Eq. (PO)) and integrating 
over the angle. Here (} = u/v and /*^°^ corresponds to the 
Fermi function (u = 0) with the chemical potential given 
by its non-equilibrium value. Similarly by substituting 
tI^^^ in dini) we obtain 



Va = 



.9/?C(TO(4 + /32) 

2w(l -/32)7/2 



(0) 
kc 



1 - f ^"^ 
^ iku 



(22) 



the momentum loss rate contribution from acoustic 
phonon scattering. 



Electronic collisions with optical phonons 
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■are 



highly inelastic. To obtain T'o, the contribution of a single 
phonon branch to the momentum loss, wc substitute 



E 



l)<5(e"'' 



k'p 



oj,)^N,5{ei;+u:,) 



(23) 

in Eq. p9)) . Here Nq ~ N{ujq) is the Bose distribution 
function evaluated at the phonon energy ujq, Wq^ is the 
golden rule expression for the transition rate from band 
a to band 7 via an interaction with a phonon of mo- 
menta q and e^^ = e^a — ^p-y- We consider the two 
optical phonon branches. Conservation of momentum 
restricts the phonon momenta to be either near the zone 
center F point or near the zone edge K point. Near 
the r point both the longitudinal and transverse opti- 
cal phonons couple to the electrons whereas near the 
K point it is mainly the A[ transverse phonon mode 
that causes inter-valley transitions. Since the typical 
phonon momentum measured from the relevant symme- 
try point is small compared to the zone boundary mo- 
menta we approximate the phonon energy by a constant: 
LUr = l96meV near the F point and ujk = 167meV near 
the K point. Furthermore we approximate w"'^ by a mo- 
menta independent constant « 2v/ {a'^y/2pLUor) for the 
zone-center phonons and g^ w 2gp for the zone boundary 
phonons21. 



B. Energy loss 

Since collisions with impurities are elastic, only 
phonons contribute to the energy loss rate 



Q = g 



2^ ^kp 



fp,)K^- 



(24) 



kap7 



The transition rate W^^ from state ka to state pj is 
given by Eg. ipg)) . 

We first consider the energy loss due to acoustic 
phonons. Electronic transitions are induced only by 
the longitudinal mode for which Wq"' = TT'D^q^{l + 
a'y cos 9) / 2 puj^. Here 9 = 6k — Op is the angle between 
the incoming and outgoing momenta, V is the deforma- 
tion potential, p is the mass density of graphene and 
LUq — cq where c is the sound velocity. We evaluate Qa 
to leading order in c/v ^ 1. Since the transitions are 
elastic to zeroth order in c/v, the leading term of Qa is 
of order of (c/v)'^. 

It is instructive to write Qa as a sum of Q™'', the 
contribution to Qa due to induced transitions, and the 
spontaneous emission contribution Qa''- The calculations 
are described in appendix 1X1 Wc find that 



Q" 



ind 



-5' 



3/32 



2pv2 (1_ ^2)7/2 



1 - f 
^ ./ku 



(25) 



As the lattice temperature increases the acoustic mode 
population increases rapidly and the rate of energy loss 
due to induced transitions increases. On the other hand, 

2^^-^^£E/^K(l-/o) + U?n (26) 
with 

/d9 
^cos(n0)/ko (27) 

depends on Tg but is independent of the lattice tempera- 
ture. At equilibrium the energy gain due to the induced 
transitions is exactly compensated by the energy loss due 
to the spontaneous emission. However, when the system 
is out of equilibrium Q^p > Q™'^ resulting in a net rate 
of energy loss by the electronic system. 

The energy loss due to the interaction of electrons with 
optical phonons Qo is evaluated using Eqs. (l23l24p . In 
expression (j23p for W^^ we make the same approxima- 
tions made above to evaluate Vo] we use a non dispersing 
phonon energy band and set w^'^ to a momenta indepen- 
dent interaction coupling constant. 

The increase in Q as the system is driven out of equilib- 
rium is due to the concomitant increase in the electronic 
temperature and in the drift velocity. However in prac- 
tice we find that since /3 <C 1 the latter contribution is 
minute so that the energy loss is dominated by heating. 

C. Numerical solution of the hydrodynamic 
equations 

Given the above expressions for the energy and mo- 
menta loss rates, we numerically solve the three cou- 
pled hydrodynamic equations. In our calculations we use 

V — 20eV for the deformation potential, and c = 0.02u 
for the sound velocity. 

1. Clean limit 

In a clean system Vi = and the momentum loss is 
only due to the phonons. In figure [2] the drift velocity is 
plotted vs. electric field for n — lO^^cm^^ and Tl — 25 
meV. Clearly du/dE is a monotonically decreasing func- 
tion of the electric field. At high fields u approaches a 
saturation value. The highest fields in this figure corre- 
spond to an electronic temperature of 200meV. 

The hydrodynamic equations imply that the drift ve- 
locity is given by Q/V. To imderstand the saturation of 
the velocity we therefore study the dependence of Q and 

V on u. In figure[3]we plot 'P{u) for different values of the 
lattice temperature. The momentum loss increases with 
E with a sharp rise at the saturation velocity Usat(21) 
indicating that the saturation of the current is due to 
the enhanced scattering at high fields. Similar behav- 
ior is obtained for the energy loss. Interestingly, partial 
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FIG. 2: Velocity saturation. The drift velocity normalized by 
the band velocity of graphene is plotted vs. electric field for 
a series of lattice temperatures expressed in meV units. 
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FIG. 3: Momentum loss rate as a function of drift velocity for 
a series of lattice temperatures expressed in meV units. 

data collapse occurs when Q and V are plotted as a func- 
tion oi fi/Te (see figure U]), demonstrating that a constant 
ratio of Q and V is reached as the neutral regime is ap- 
proached. 

We find that u^at is only weakly density dependent 
in the range n = 0.1 — 10 • lO^^cm"^. Therefore the 
saturation current Isat = enugat is, to a good approx- 
imation, linear in n. We also find that for all values 
of the applied field Pq ^ Pa- Acoustic phonons there- 
fore play an essential role in the non linear dc electrical 
properties of graphene. The energy loss rate in doped 
graphene is dominated by acoustic phonons at low elec- 
tronic temperatures-^^ , however far from equilibrium it is 
the optical phonons that dominate Q. 

In deriving Vo and Qo we implicitly assumed that the 
optical phonons are thermalized. Experimental and the- 
oretical work has shown that the optical phonons may be 
far from equilibrium in hot carbon nanotubeSj^i^ and 
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FIG. 4: The momentum loss rate as a function of /^/Te- Dif- 
ferent curves correspond to different values of the lattice tem- 
perature in meV units. Data collapses as the neutral regime 
is approached. The drift velocity is limited by the ratio of the 
typical energy loss to the typical momentum loss averaged 
over transitions, which is proportional to the Dirac band ve- 
locity. 

this is also a possibility in graphene. The approximation 
we make here in setting the temperature of the optical 
phonons to the temperature of the acoustic phonon bath 
Tl can be justified a posteriori by our numerical results. 
We find that Q^'° ^ Q™'^ far from equilibrium. Any in- 
crease in the temperature of the optical phonon bath will 
infiuence only Q*"'* and will thus have little effect on our 
non-linear transport results. 

D. Influence of disorder 

The presence of disorder leaves the form of the hydro- 
dynamic equations unchanged, however it does increases 
the momentum loss rate V. Therefore, for a given value 
of the electric field, we expect the drift velocity u = Q/V 
to be reduced relative to its value in a clean system. 
In figure O the drift velocity is plotted as a function of 
the electric field for a graphene sheet with an electronic 
density n = 10^^cm~^ and an impurity concentration 
Hi = lO^^cTO"^. For the strongest fields = 280meV. 
Comparing figure [5] and figure [2] clearly shows that the 
drift velocity is indeed reduced. Moreover we find no cur- 
rent saturation even at the highest fields. Since disorder 
does not modify the energy loss, we expect saturation to 
occur only when Pi. 

V. DISCUSSION 

At high electronic temperatures, rapid e-e collisions 
justify the hydrodynamic theory used here. The high 
temperatures are characteristic of the non-linear response 
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FIG. 5: Current as a function of electric field for n = 
10^^ cm~^ and rii = 10^^ cm~^. The different curves corre- 
spond to difi'erent values of the lattice temperature in meV 
units. 

regime which is our principle focus, but the hydrody- 
namic theory will also apply in the linear regime for suf- 
ficiently clean and hot samples. 

In the linear response regime we find that inter-particle 
collisions have little effect on the resistivity when fi/T is 
large. However at high temperatures, or alternatively at 
low densities, interactions between electrons increase the 
resistivity. Recently Bolotin et. al. measured the re- 
sistivity in ultra-clean suspended graphene over a wide 
range of densities and temperatures^'*. Their experimen- 
tal samples were close to the ballistic limit in which the 
quantum contact resistances are a substantial fraction of 
the overall resistance. Surprisingly this study found that 
the phonon induced resistivity is density dependent at 
low densities. Although our theory can not make quanti- 
tative predictions in the ballistic regime since it is unable 
to account for the quantum resistance, our results do sug- 
gest e-e interactions as a possible origin of this density 
dependence. 

In the non-linear regime we find that u = Q/V: the 
drift velocity is given by the ratio between the energy 
loss rate and the momentum loss rate. At strong fields 
high electronic temperatures are responsible for a rapid 
increase of both Q and V that results in the saturation 
of M at a velocity of the order of 10^ cm / sec. The satu- 
ration velocity is only weakly density dependent there- 
fore the saturation current I^at — enusat mcreases, to 
a good approximation, linearly with density. The elec- 
tronic temperature at which velocity saturation occurs 
increases with the impurity concentration. For even mod- 
erately disordered samples, electronic temperatures reach 
unphysically large values before saturation occurs. 

In this work we neglected the phonons of the substrate. 
The interaction of Dirac quasi-particles with substrate 
phonons can however be important, depending on the 
type of substrate and its distance from the graphene 
sheet^i^. Following the prescription given in section 
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IIVI the contribution of these phonons may straightfor- 
wardly be added to the the hydrodynamic equations 
(|5|7|) . Strong couphng between the substrate phonon 
modes and the electrons may significantly lower the 
electronic temperature in which current saturation oc- 
curs. If the coupling to substrate phonons is too strong 
{Tph < Tee), the hydrodynamic approach is invalidated. 

Previous work attributed current saturation in carbon 
nanotubes^ to the the sudden onset of momentum re- 
laxation by zone-boundary optical phonons. More re- 
cently Meric et. al. associated the saturation of the cur- 
rent in graphene based field-effect transistors with surface 
phonons of the Si02 substrate^. In this work we find 
an additional mechanism for current saturation in which 
the increase of the drift velocity is inhibited by high elec- 
tronic temperature at strong fields. The two mechanisms 
can be distinguished experimentally by the dependence of 
saturation velocity on carrier density which they predict. 
The electron-heating mechanism leads to a saturation ve- 
locity which is ~ 10% of the Dirac velocity and weakly 
dependent on carrier density. The saturation current is 
therefore, to a good approximation, proportional to car- 
rier density. The phonon back-scattering mechanism, on 
the other hand, predicts a saturation velocity which is 
~ VLUph/n and therefore a critical current which varies 
as the square root of carrier density. 

In real devices the top-gate lies in close proximity 
to the source-to-drain conduction channel. The electric 
field in the channel is consequently strongly space depen- 
dent. The hydrodynamic theory outlined above may be 
straightforwardly generalized to include such spatial vari- 
ations. The out of equilibrium system is then described 
by three hydrodynamic functions u(r), Te{r) and /i(r) 
whose values are determined by the continuity equation, 
the momentum balance equation and the energy balance 
equation. 

After completion of this paper we became aware of an 
interesting and closely related recent study by Barreiro 
et. al.i^ This complementary work addresses high field 
transport in graphene both experimentally and theoreti- 
cally. The experimental results support our findings re- 
garding the absence of current saturation in disordered 
graphene. The theoretical analysis in Ref.(27) utilizes 
the Boltzmann equation, however it neglects e-e inter- 
actions and the interactions of the Dirac quasi-particles 
with acoustic phonons. As we explained above e-e in- 
teractions are expected to be important when the elec- 
tronic temperature is high, and interactions with acoustic 
phonons are essential for the description of current satu- 



ration in clean graphene. 
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APPENDIX A: ENERGY LOSS TO ACOUSTIC 
PHONONS 

The energy loss rate due to the interactions of the elec- 
trons with acoustic phonons 



Qa 



2pc 



k,pa7 

(N, + l)5{et; - u;,) + N,6{et; + u;,)] (Al) 



where q — |k — p|, is evaluated to leading order in 
f3 — c/v. Such an expansion is valid in the quasi-elastic 
scattering regime when Te > Tbg- 

The energy conservation condition in (jAl|) inhibits 
inter-band transitions. Furthermore it sets the value of 
the dummy variable p to 



Po^k± 2akP\ sin(6i/2)| + 2/3^ sm'^{0/2) 



(A2) 



where the positive sign relates to a phonon emission pro- 
cess and the minus sign to a phonon absorption process. 
Hence to first order in f3 



S{et;±c\k-p\ 



^a /^^ [lT«7hin(g/2)|]. 



(A3) 

The total energy loss (jAip is a sum of Q^p the energy 
loss due to spontaneous emission of phonons and Q™'' 
the energy gain due to induced transitions. For the latter 
Eqs. l|AllA2IA3|) imply 



Qtnd ^ 



8V^T,. f k^dkd9kd0. 



(27r)3 



I 

^sinHY^^hc- (A4) 



Integrating first over Op and then over 0k results in 
expression (|25|l for Q™''. A similar derivation yields 
Eq.®. 
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